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1.  INTRODUCTION 

This  annual  report  covers  research  concerned  with 
increasing  our  understanding  of  nonideal  dense  plasmas 
produced  by  strong  explosively  driven  shock  waves.  The 
research  effort  was  performed  for  the  Office  of  Naval 
Research-Power  Program  under  Contract  N00001i<-78-C-n35l| , 
and  utilized  a novel,  high  performance  explosively  driven 
nlasma  source. 

The  Dlasma  conditions  produced  by  these  new  sources 
are  substantially  more  energetic  and  nonideal  than  those 
investigated  or  achieved  in  earlier  programs,  and  under- 
score the  dearth  of  knowledge  and  understanding  of  dense 
nonideal  plasmas. 

The  research  effort  covered  in  this  report  consists 
mainlv  of  four  explosively  driven  plasma  source  experiments 
that  were  designed,  fabricated,  instrumented,  fired,  and 
analyzed  during  the  course  of  this  contract  period. 

1.1  Background 

\ 

For  several  years  Artec  Associates  has  been  actively 

9 

enraged  in  pulsed  power  generation  using  an  explosively 
driven  MHD  concept-1-  * 2 ’ 3 . The  key  element  of  this  novel 
concept  is  a nlasma  source  which  converts  the  chemical 


energy  of  an  explosive  into  a very  high  enthalpy  fully 
ionized  plasma.  The  high  density  plasma  is  in  turn  driven 
into  a channel  and  through  a magnetic  field  in  a MHD  section, 
'’’he  electrical  output  of  this  section  is  used  to  power  an 
augmenting  magnetic  field  intensifying  the  total  magnetic 
field  through  self-excitation  and  ultimately  driving  an 
external  load. 

The  concept  of  a fully  self-excited  explosive  MHD 
pulse  nower  eenerator  is  particularly  attractive  for  many 
Naval  applications . Such  a venerator,  consisting  only  of 
explosive,  metal  tubing  and  a permanent  magnet,  would  have 
an  attractively  long  shelf-life.  Startup  time  from  initia- 
tion to  full  nower  is  nearly  instantaneous,  and  tens  of 
mega.loules  of  electrical  energy  are  available  at  gigawatt 
n'lw^r  levels.  These  explosive  generators  are  Inherently 
small,  and  conoatible  with  many  special  weaoon  applications. 
Other  applications  may  well  become  apparent  as  the  program 
progresses . 

Fxnlosive  MHD  generators  of  the  externally-excited 
varietv  have  been  experimentally  demonstrated  many  times 
in  the  past,  but  the  real  progress  toward  a fully  self- 
excited  generator  has  been  thwarted  by  a lack  of  knowledge 
of  the  properties  of  dense  nonideal  plasmas.  Measured 
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electrical  conductivities  In  our  experiments  for  example  are 
substantially  less  (up  to  a factor  of  10)  than  calculations 
based  on  Ideal  plasma  theory.  To  compensate  for  the  difference 
and  to  still  achieve  full  self-excitation  it  has  been  necessary 
to  increase  both  enemy  and  duration  of  the  explosive  plasma 
source. 

1.2  Previous  Research  Effort 

Previous  research  has  been  directed  toward  investigating 
the  s^ock  nhvslcs  of  nonideal  plasmas^-,  and  developing  an 
advanced  hirh-perfornance  nlasma  source.  This  effort  was 
ouite  successful,  and  Artec  Associates  conceived  an  exolo- 
- i ve  nlasma  source  that  was  far  more  energetic  than  anv 
used  before. 

The  new,  hi  oh  performance  explosive  plasma  source 
connr^sses  plasma  to  denser,  no^e  nonideal  states  than 
before,  'duch  of  the  enerfv  addition  to  the  plasma  occurs 
at  densities  exceeding  normal  solid  densitv  for  many  mate- 
rials and  nonideal  effects  dominate  the  plasma  behavior. 
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2.  EXPERIMENTAL  RFSFARCH  ON  NON I DEAL  PLASMAS 
2.1  Descriotion  of  the  Plasma  Source 

A schematic  diagram  of  the  plasma  source  Is  shown  In 
Figure  1.  It  is  comprised  of  an  annular  region  of  argon 
gas  which  is  shock-compressed  by  a detonating  explosive. 

At  the  end  of  the  cylindrically  symmetric  annular  region, 
the  shocked  argon  is  converged  radially  into  a hemispherical 
transition  section  and  exhausted  into  a diagnostic  channel 
located  on-axis.  This  driver  configuration  was  used  in 
three  shots,  130-2,  130-3>  and  130-4.  The  conical  driver 
shown  in  Figure  2 and  discussed  in  previous  reports  was 
used  to  drive  Shot  130-1. 

The  annular  portion  of  the  plasma  source  operates 
much  as  the  explosive  plasma  sources  described  in  earlier 
research  reoorts1  > 2 » 3.  More  energy  is  add-ed  to  the  plasma 
by  the  pronounced  radial  convergence  of  the  flow  and  explo- 
sive compression  of  the  reservoir.  In  other  words,  the 
converging  shock  flow  now  stagnates  along  the  axis,  forming 
a much  higher  pressure,  higher  density  plasma  reservoir  than 
previously.  Nearly  simultaneously,  the  metal  caD  at  the  end 
of  the  interior  explosive  is  explosively  accelerated  and  com- 
presses and  energizes  the  plasma  even  further.  As  can  be 
annreciated,  the  convergence,  compression,  and  exhaust 
processes  are  all  mutually  dependent.  Timing  is  critical 
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in  achieving  high  performance  of  the  plasma  source  but 
experimental  performance  of  the  plasma  source  has  been 
very  satisfying,  even  though  the  design  has  not  yet  been 
optimized . 

Overall  efficiency,  i.e.  conversion  of  explosive  chem- 
ical energy  to  useful  plasma  enthalpy,  already  exceeds  25$, 
and  will  most  likely  be  increased  further  with  design  re- 
finements. We  note  that  plasma  energy  is  higher  than  that 
achieved  in  any  previous  plasma  source.  MHD  channel  con- 
ditions are  also  very  good,  and  even  with  a relatively 
small  diameter  driver,  a significant  degree  of  generator 
self-excitation  has  been  achieved.  Full  self-excitation 
appears  highly  probable  with  a larger  source,  or  with 
mechanical  design  refinements. 

2.2  Estimate  of  Monideal  Plasma  States 

Plasma  states  achieved  in  the  new  plasma  source  may 

be  estimated  from  experimental  data  and  the  best  available 

4 

nonideal  plasma  theory  . The  initial  driver  shock  wave  in 
the  annular  region  of  the  plasma  source  has  a measured  par- 
ticle velocity  of  8.5  km/s  (the  detonation  velocity  of 
Octol).  Using  this  measured  velocity  and  assuming  the 
Debye-Huckel  theory  with  a zero  shielding  coefficient,  the 
shocked  gas  conditions  are  calculated  to  be:  pressure} 
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3GPa  ( 30  kbar);  density,  0.33  Mg/m^  (g/cm^);  and  energy 
37  MJ/kg.  These  estimates  are  fairly  accurate,  as  the 
pressure  and  energy  states  behind  a shock  are  relatively 
Insensitive  to  uncertainties  in  the  equation  of  state. 
However,  the  Debye-Huckel  plasma  theory  used  in  the  cal- 
culations is  about  at  the  limit  of  credibility  because  of 
the  dense  plasma  state. 

An  estimate  of  the  reservoir  entropy  state  may  be 
calculated  assuming  a single  strong  reflected  shock.  To 
account  for  the  high  density,  we  have  used  the  empirical 
Jones-VJilklns-Lee  (JWL)^*^  equation  of  state  commonly  used 
in  explosive  detonation  theory.  Calculated  conditions  for 
the  above  initial  shock  are:  pressure,  35  OPa  (350  kbar); 
density,  1.30  Mg/m^  (g/cm^);  and  energy,  80  MJ/kg. 


Final  compressed  state  in  the  reservoir  can  be  estimated 
assuming  an  isentropic  compression  of  this  stagnated  region 
to  a value  obtained  by  matching  experimental  shock  velocity 
data.  For  example,  the  magnitude  of  the  shock  velocity  in 
the  diagnostic  channel  is  a direct  consequence  of  specific 
reservoir  conditions.  Rv  matching  the  experimental  data 
for  shock  velocity,  a unique  value  of  the  reservoir  condi- 
tions can  be  obtained.  Using  this  approach,  peak  plasma 
conditions  in  the  reservoir  were  determined  to  be:  pressure, 
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8^  OPa  (8*10  kbar):  density,  2.2  Mg/m^  (g/cm3);  and  energy 
97  MJ/kg. 

Plasma  states  in  the  reservoir  achieved  by  the  new 
plasma  source  are  located  almost  off  the  edge  of  the  plasma 
parameter  chart,  shown  in  Figure  3,  and  even  the  expanded 
plasma  is  well  into  the  nonideal  region  at  the  right  of 
the  plasma  parameter  chart.  The  chart  and  the  parameters 
were  described  previously  in  our  1978  research  report^-, 
but  key  parts  will  be  summarized  here.  QNTM  is  the  ratio 
of  electron  thermal  deBroglie  wavelength  to  Debye  length. 

For  values  greater  than  0.1  quantum  diffraction  effects 
become  important  in  plasma  theory.  INTPAR  is  the  plasma 
nonideal  interaction  parameter,  i.e.  the  ratio  of  coulomb 
electrostatic  energy  to  thermal  kinetic  energy.  Ideal 
plasma  theory  assumes  INTPAR  to  be  very  snail.  In  fact 
the  Debye-Huckel  theory,  which  represents  a first  order 
correction  to  ideal  plasma  theory,  is  known  to  be  invalid 
above  an  interaction  parameter  of  0.^7.  Finally,  DON  is 
the  quantum  degeneracy  parameter,  the  ratio  between  electron 
thermal  deBroglie  wavelength  to  electron  separation. 

It  can  be  seen  from  the  figure  that  the  Artec  plasma 
source  yields  conditions  well  into  the  nonideal  plasma 
regime,  indeed  much  further  than  the  earlier  plasma  sources. 
In  order  to  understand  the  details  of  the  shock  energization 
process,  particularly  the  final  compression  of  the  reservoir 

-17- 
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Fig.  3 Plasma  Parameter  Chart  for  Argon 
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it  is  necessary  to  have  a clear  understanding  of  the  physics 
of  strongly  nonideal  plasmas.  Once  the  nlasma  expands  into 
the  MTO  channel,  the  static  pressures  drop  to  about  0.1  GPa 
(1  kbar),  hence  the  MHD  interaction  process  can  be  studied 
using  "moderately  nonideal"  plasma  theory. 

2.3  Description  of  the  Experimental  Apparatus 

Pour  experiments  were  designed,  instrumented  and  fielded 
for  this  contract.  The  basic  design  parameters  are  Riven  in 
Table  1. 

Shot  130-1  utilized  a conical  driver  and  xenon  gas  at 
137  kilopascals  (20  psla)  as  the  driver  gas.  The  main 
ournose  or  this  experiment  was  to  fire  a conical  imploding 
driver,  measure  the  flow  properties  of  the  driver  gas  as  it 
expanded  into  a partiallv  evacuated  channel,  and  compare 
these  properties  with  the  expected  flow  properties  based 
on  calculated  driver  performance. 

Shot  130-2  utilized  a 17  cm  diameter  imploding/exploding 
driver  with  xenon  gas  at  6?^  kilopascals  (100  psia)  as  the 
driver  gas.  The  main  purposes  of  this  experiment  were  to 
evaluate  this  type  of  driver,  to  measure  the  properties  of 
the  plasma  as  it  expanded  into  a partially  evacuated  channel, 
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and  to  assess  the  improvement  in  performance  of  this  source 
by  comparing  the  output  with  the  configuration  of  Shot  130-1. 

Shot  130-3  utilized  a 17  cm. diameter  imploding/exploding 
driver  as  before  but  with  argon  at  2.12  megapascals  as  the 
driver  gas.  The  purpose  of  this  experiment  was  to  fire  an 
imploding/exploding  type  plasma  source  using  argon  and  com- 
pare the  resulting  performance  with  the  same  configuration 
using  xenon  (Shot  130-2).  Plasma  output  was  to  be  assessed 
using  the  same  diagnostics  as  Phot  130-2,  and  in  addition 
a self-excited  MHD  station  was  included  for  additional  diag- 
nosis of  the  plasma  parameters. 

Shot  130-4  utilized  a 30  cm.  imploding/expioding  driver 
with  argon  at  2.07  megapascals  as  the  driver  gas.  It  was 
relt  that  the  mass  of  the  driver-processed-gas  in  the 
17  cm.  diameter  drivers  was  limiting  the  flow  duration  of 
the  channel  flow  to  nominally  30  microseconds,  and  that 
increasing  the  outside  diameter  to  30  centimeters  and  the 
mass  of  the  gas  by  better  than  fivefold  would  increase  the 
total  flow  time  to  50  to  60  microseconds,  and  permit  as  much 
as  40  dB  magnetic  field  enhancement. 

Additionally,  a minor  redesign  of  the  field  coil  to 
prevent  mechanical  short  circuits  was  also  included  on  130-4 
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to  negate  suspected  shorts  which  we  initially  suspected 
heloed  to  limit  th®  increase  in  magnetic  field  in  130-3 
to  iQ  dn. 


2. 4 Description  of  the  Diagnostics 

Shot  130-1:  A drawing  of  the  Dlasna  source  used  in 
this  shot  was  shown  previously  in  Fi cure  2.  The  explosive 
used  was  75/25  Octol;  the  Imploding  liner  was  copper  (4.0  mm. 
cr  156  mils  thick);  the  cone  half-angle  was  10°;  and  the 
charge-to-mass  ratio  was  3.0.  The  thick  steel  wave  shaper 
permitted  a 3ingle  detonator  to  provide  a symmetrical  annular 
detonation  front.  The  cone  angle  was  calculated  to  produce 
an  axiallv  progressing  impact  point  with  a speed  of  17  km/s. 


The  driver  gas  chamber  was  loaded  with  20  grams  of 
xenon  to  provide  a nominal  shock  pressure  of  2.4  GPa  (24 
kbar).  The  diagnostic  channel  was  loaded  with  air  at  6.65 
kPa  (50  Torr)  and  was  separated  from  the  driver  by  a thin 
steel  diaphragm. 

The  layout  of  the  diagnostics  in  the  channel  is  shown 
in  Figure  4.  There  was  one  velocity  gage  or  open  circuit 
voltage  measurement  and  one  electrode  station  with  an  11 
milliohm  load  (a  thin  strip  of  stainless  steel)  at  which 
the  loaded  voltage  and  currer.t  were  to  be  measured.  A pair 
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of  0.3  meter-long  barium- ferrite  magnets  provided  a trans- 
verse magnetic  field  of  77  mT  (770  gauss)  and  33  mT  (830 
gauss)  at  the  velocity  gage  and  load  stations  respectively. 
There  were  two  conductivity  papres,  one  located  near  the 
upstream  end  of  the  magnets  and  one  near  the  downstream 
end  as  shown  in  the  figure. 

Shot  130-2:  A drawing  of  the  plasma  source  was  shown 
previously  in  Figure  1.  The  explosive  used  was  75/25  Octol. 
The  wave  shaper  was  so  configured  that  a single  detonator 
on-axis  provided  a symmetric  detonation  front  sweeping 
along  both  tubes  simultaneously. 

The  driver  gas  chamber  was  loaded  with  56  grams  of 
xenon  to  provide  a nominal  shock  pressure  of  3*0  GPa  (30 
kbar).  The  diagnostic  channel  was  loaded  with  air  at  6.65 
kPa  (50  Torr)  and  was  separated  from  the  driver  by  a 25pm 
thick  stainless  steel  diaphragm. 

The  layout  of  the  diagnostics  in  the  channel  is  shown 
in  figure  4 and  is  identical  with  Shot  130-1.  There  was 
one  velocity  gage  or  open  circuit  voltage  measurement  and 
one  electrode  station  with  an  11  milliohm  load  (a  thin 
strip  of  stainless  steel)  at  which  the  loaded  voltage  and 
current  were  measured.  A pair  of  0.3  meter-long  barium- 
ferrite  magnets  provided  a transverse  magnetic  field  of 
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75  mT  (7^9  gauss)  and  89  rnT  (890  cjauss)  at  the  velocity 
gage  and  load  stations  respectively.  There  were  two  con- 
ductivity gashes,  one  located  near  the  upstream  end  of  the 
magnets  and  one  near  the  downstream  end  a3  shown  in  the 
figure. 
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Figure  5 Layout  of  Self- Excited  Generator  Station 
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coll  were  35.9  nanohenries  and  0.^03  mllllohms  respectively. 
Measured  D/I  ratio  at  the  tube  wall  was  18.2  mlcroTesla/amp . 

Shot  130-JI:  The  geometrical  configuration  of  130-^ 
was  nearly  Identical  to  130-3  except  that  the  driver  outside 
diameter  was  Increased  to  30  centimeters.  A full  2X  scale 
up  from  the  baseline  17  centimeters  was  discarded  since 
tubing  of  the  desired  dimensions  is  not  commercially  avail- 
able. The  driver  gas  was  argon  and  the  diagnostic  channel 
was  loaded  with  air  at  6.^5  kiiopascals  (50  Torr)  as  pre- 
viously . 

The  first  section  of  the  diagnostic  channel  contained 
2 side-mounted  conductivity  gages  (gages  1 and  2),  a 
velocity  gage  and  a conventional  MHD  generator.  The  mag- 
netic field  for  these  3 diagnostic  tools  was  provided  by 
two  barium-ferrite  permanent  magnets  mounted  to  produce  a 
D field  of  .09  Tesla  orthogonal  to  the  flow  velocity  vector. 

The  second  section  of  the  diagnostic  channel  included 
a third  side-mounted  conductivity  gage  and  a self-excited 
MKP  generator  stage.  The  experimental  configuration  is 
given  schematically  In  Figure  6. 

The  Initial  magnetic  field  for  this  section  was  pro- 
vided by  a pair  of  barium-ferrite  permanent  magnets. 
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Additionally,  flush  brass  electrodes  were  connected  to  a low 
inductance  one-turn  self-excitation  coil  made  from  # 8 AWG 
wire  braid.  The  measured  inductance  of  the  field  coil  was 
37.9  nanohenries.  The  D.C.  resistance  was  0.525  milliohms. 
The  measured  B/I  at  the  tube  centerline  was  24.3  microtesla 
per  amp. 
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3.  RESULTS  AND  CONCLUSIONS 
3.1  Plasma  Source  Performance- 

Pairs  of  symmetrical  ionization  pins  placed  on  either 
side  of  the  plasma  source  reported  within  fractions  of  a 
microsecond  of  each  other,  verifying  the  symmetry  of  the 
detonation.  The  pins  also  ”erified  the  correct  (expected) 
detonation  velocity  for  Octol,  8.48  km/sec.  A typical  plot 
of  driver  ionization  pin  data  is  shown  in  Figure  7-  The 
precision  and  symmetry  of  the  data  is  clear  from  the  figure. 

In  all  cases  the  driver  p:as  initial  pressure  (hence 
density)  was  adjusted  to  yield  30  kbars  static  pressure 
behind  the  driven  initial  shock.  This  nominal  or  standard 
condition  was  used  in  each  case  to  allow  comparing  the 
different  cases,  one  with  another. 

A flash  X-ray  photograph  of  an  annular  plasma  source, 
shown  with  the  detonation  wave  approximately  half  way  along 
the  driver,  is  given  in  Figure  8.  The  high  degree  of  sym- 
metry both  for  the  inside  and  outside  collapsing  walls  and 
from  one  side  to  the  other  is  evident. 

In  the  linear  section  of  the  source  for  130-2  and  130-3j 
the  total  enthalpy  was  80.2  MJ/kg  giving  a total  energy  of 
4.5  mj  for  the  86  grams  of  gas.  In  the  convergent  section 
of  the  source,  which  effectively  acts  as  a reservoir,  the 


Distance  - Time  Diagram  for  Octol  Annulor  Plasma  Source  (130-4) 


total  enthalpy  Is  estirated  to  be  13*1  MJ/kg.  On  the  assump- 
tion that  all  of  the  gas  achieved  this  enthalpy,  the  total 
energy  In  the  plasma  was  7.5  MJ. 

3.?  Channel  Conditions 

The  initial  pressure  in  the  diagnostic  channel  in  all 
cases  waa  6.65  kilopascals  (50  Torr)  of  air,  hence  the 
performance  of  all  the  sources  in  terms  of  the  equivalent 
reservoir  conditions  Droduced  can  be  seen  from  a comparison 
of  the  initial  shock  velocitv  in  the  diagnostic  channel. 

^hese  velocities  are  tabulated  in  Table  2. 

^or  Shot  130-1,  a channel  shock  speed  of  27-3  km/s 
was  obtained  from  shock  front  arrival  times  at  the  diagnostic 
stations  and  from  shock  breakout  at  the  diaphragm  which  was 
monitored  by  a photodiode.  The  initial  plasma  flow  velocity 
calculated  from  the  observed  shock  velocity  is  24.8  km/s. 

Prom  the  velocity  gage  record  the  initial  flow  velocity 
was  measured  to  be  22.7  km/s,  which  is  S percent  below  the 
calculated  value.  Since  this  particular  velocity  gage  was 
located  five  channel  diameters  downstream  of  the  magnet  edge 
and  3ince  the  magnetic  Reynolds  number  of  the  flow  was  not 
sufficiently  high  to  cause  significant  perturbation  of  the 
velocity  gage  reading  at  that  distance,  we  suspect  the  source 
of  the  discrepancy  must  lie  in  other  than  Rm  effects. 


TABLE  2 


INITIAL  SHOCK  AND  FLOW  VELOCITIES 


Measured  Initial 
Shock  Velocity 
Shot  No . ( km/sec ) 


Measured  Initial 
Flow  Velocity 
( km/sec ) 


Calculated  Initial 
Flow  Velocity 
(km/sec ) 


130-1 

27.3 

22.7 

24.8 

130-2 

25.9 

19-2 

23-6 

130-3 

26.4 

20.6 

23.6 

130-4 

32.5 

270 

29.7 
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For  Shot  130-2,  a channel  shock  speed  of  25.9  km/s 
was  obtained  from  shock  front  arrival  times  at  the  diagnos- 
tic stations  and  from  shock  breakout  at  the  diaphragm,  which 
was  monitored  by  a photodiode.  This  slightly  lower  velocity 
(5  percent  compared  to  the  nrevious  shot)  is  probably  within 
the  experimental  scatter  in  the  data. 

For  Shot  130-3»  a channel  shock  velocity  of  26.^  km/s 
was  measured,  compared  to  25.9  km/s  in  the  previous  shot, 
and  is  essentially  the  same,  implying  that,  from  a shock 
velocity  point  of  view,  argon  and  xenon  drivers  are  nearlv 
identical  in  performance. 

A channel  shock  velocity  of  32.5  km/s  was  measured 
for  Shot  130-^  compared  to  2 6.4  kn/s  and  25.9  km/s  for 
the  previous  shots,  respectively.  Hence  it  is  clear  that 
the  larger  driver  and  its  increased  convergence  of  driver 
gas  oroduce  greater  peak  reservoir  conditions,  hence  greater 
shock  velocities  at  the  diagnostic  stations.  Indeed,  this 
shock  velocity  is  the  highest  that  has  been  obtained  in 
the  program  to  data,  and  as  such  represents  a milestone  in 
our  progress. 

Tabulated  in  Table  2 are  the  measured  and  calculated 
initial  f]ow  velocities  for  these  four  experiments.  One 


Flow  Velocily  (km/s) 
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velocity  histories,  both  in  absolute  value  of  velocity  and 
flow  tine. 

The  30  cm.  diameter  driver  of  130-4  reives  substantially 
greater  flow  velocities;  indeed  the  highest  peak  obtained 
in  our  research  to  date.  These  velocities  are  greater  at 
any  given  time  and  are  maintained  to  approximately  40-45 
microseconds  in  time.  The  efficacy  of  the  larger  diameter 
driver  is  clearly  apparent  from  this  trace.  However,  the 
decay  of  the  channel  is  greater  than  expected.  The  decline 
in  velocity  implies  that  the  higher  reservoir  conditions  are 
not  being  maintained  as  long  as  would  have  been  expected 
with  the  additional  mass  of  working  gas  produced  by  the 
larger  driver. 


These  points  suggest  the  possibility  that  other  pre- 
viously unsuspected  phenomena  may  be  limiting  te^t  flow 
time  both  for  the  small  driver  as  well  as  the  large.  A 
prime  candidate  is  the  premature  "pinching  off"  of  the 
flow  at  the  channel  entrance  by  distortions  of  the  entrance 
region  by  the  very  high  reservoir  pressures.  This  type  of 
behavior  has  been  observed  by  others,  for  example,  the  com- 
puter experiments  of  the  Voitenko  compressor  performed  by 

O 

Brown  and  Lohmann  of  LLL  . They  observed  that  at  very 

high  reservoir  pressures,  well  in  excess  of  the  yield  stresses 
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In  the  reservoir  wall  material,  the  reservoir  end  wall  near 
the  channel  is  extruded  in  the  channel  entrance  region, 
"pinching  off''  the  outflow.  Even  if  the  flow  is  not  com- 
pletely blocked,  a "nozzle"  is  created  at  the  channel 
entrance  requiring  the  gases  expand  further  in  a steady 
expansion  before  shifting  to  an  unsteady  expansion.  As 
the  steady  expansion  is  less  efficient  in  translating 
pressure  into  velocity  than  the  unsteady  expansion  for 
flow  Mach  numbers  greater  than  unity,  less  shock  velocity 
is  obtained  for  a given  reservoir  condition.  A discussion 
of  this  subtle  but  not  commonly  appreciated  point  is 
given  in  Reference  7. 

This  problem,  and  modlf ications  to  the  reservoir — 
channel  transition  area  to  counter  it,  will  be  addressed 
in  future  work. 

3.3  Expanded  Plasma  Electrical  Properties 

The  electrical  conductivity  of  the  gas  in  the  channel 
is  of  key  interest  in  this  research.  Not  only  are  precise 
values  of  conductivity  presently  impossible  to  obtain  theo- 
retically due  to  the  complex  nonideal  nature  of  the  plasmas, 
but  the  measurements  also  include  significant  scatter  in 
results,  both  from  similar  instrumentation  and  between 
different  measurement  techniques. 
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For  example,  Figure  10  shows  electrical  conductivity 
determined  from  an  eddy  current  gage  (described  earlier)  and 
that  determined  from  examining  the  V-I  characteristics  of 
the  loaded  generator  section.  One  notes  that  the  conduc- 
tivities vary  by  about  a factor  of  two  with  the  eddy  current 
gage  giving  the  higher  values. 


The  conductivity  results  from  Shot  130-2  are  shown  in 
Figure  11.  In  this  case,  two  electrical  conductivity 
gages  were  used,  one  upstream  and  one  downstream  of  the 
load  section  (see  Figure  4).  The  gages  gave  similar  results 
with  the  downstream  gage  giving  generally  slightly  higher 
values.  One  also  notes  that  the  conductivity  histories  are 
maintained  in  time  better  than  130-1  as  would  be  expected 
from  examining  the  particle  velocity  histories  of  the  two 
shots . 

As  previously,  one  also  notes  that  the  values  determined 
from  the  load  section  are  lower  by  about  50  percent. 


The  electrical  conductivity  data  for  Shot  130-3  are 
given  in  Figure  12.  In  this  shot,  3 eddy  current  gages 
were  used  to  obtained  conductivity  in  addition  to  the 
valuer  determined  from  the  load  resistance,  and.  self-excited 
generator  sections.  The  locations  of  the  5 different 
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measurements  can  be  determined  from  Figure  6 and  should  not 
be  dismissed  in  accounting  for  the  different  conductivity 
values . The  time  values  shown  are  from  the  start  of  each 
trace  and  are  at  different  absolute  times.  A plot  of  the 
values  reordered  to  a laboratory  time  scale  shows  some 
interesting  effects.  For  example,  in  Figure  13  we  note 
that  conductivity  appears  generally  to  be  increasing  as 
the  flow  progresses  downstream  and  has  a generally  decreas- 
ing trend  with  time  at  a given  measurement  station. 

We  also  note  that  the  conductivity  measurements  based 
on  the  load  measurements  and  self-excited  generator  measure- 
ments show  an  increase  in  conductivity  for  the  first  few 
(<  5)  microseconds.  If  this  is  a transient  effect  caused 
by  say,  electrode/gas  processes,  then  the  agreement  in  the 
shape  begins  to  look  better.  The  last  maxima  in  the 
generator  derived  trace  remains  unexplained.  As  the  traces 
for  130-4  are  unreliable  late  in  the  flow,  no  additional 
data  is  available  for  comparison. 

We  expect  that  more  might  be  gleaned  from  this  data 
by  an  analysis  of  the  various  data  channels  microsecond 
by  microsecond,  and  propose  to  use  such  a technique  on 


future  runs . 
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Figure  14  Conductivity  Comparison  for  Xenon  and  Argon  Shots  (130-2,  130-3) 
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Figure  15  Electrical  Conductivity  Histories,  Shot  130-4 
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For  example,  field  coll  data  at  the  generator  station 
for  130-3  are  shown  in  Figure  16.  Belf-excitation  is  occur- 
ring throughout  the  entire  pulse,  as  evidenced  by  increasing 
current  and  positive  voltage.  A short  circuit  occurs  at 
3^us,  probably  due  to  failure  of  the  Mvlar  plastic  insula- 
tion in  the  coil.  Total  magnetic  field  amplification  during 
the  pulse  was  10  dB  (a  factor  of  9).  Initial  increase  was 
1.3  dB/ps,  falling  to  0 . h dB/ps  at  the  end  of  the  pulse. 

The  MHD  generator  was  analyzed  for  consistency  by 
using  measured  plasma  parameters  to  calculate  performance, 
and  then  comparing  them  with  observed  values.  The  B-field 
magnitude  was  calculated  from  measured  initial  field,  current, 
and  B/I . This  was  combined  with  velocity  to  yield  the  open 
circuit  Faraday  voltage  (Bub)  shown  in  Figure  17. 

Initially,  at  I = 0,  the  calculated  open  circuit  voltage 
exceeds  measured  field  coil  voltage.  This  may  be  caused  by 
internal  inductance  of  the  generator,  forming  an  inductive 
voltage  divider  with  the  load.  A value  of  16  nanohenries 
agrees  with  the  data,  but  is  larger  than  theoretical  esti- 
mates, which  yield  about  5 nanohenries.  On  the  other  hand, 
internal  inductance  of  the  load  station  had  previously 
been  determined  to  be  14  nanhenries.  An  alternate  explana- 
tion is  loading  of  the  generator  by  eddy  currents,  expected 
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because  of  the  very  short  coil  length;  this  loading  would 
act  much  like  an  inductive  divider  in  the  present  experiment . 
Some  reduction  of  open  circuit  Bub  was  expected,  but  the 
observed  amount  of  about  3355  appears  excessive.  Better 
design  of  the  field  coil  assembly  might  well  be  able  to 
reduce  it  to  15$  or  less. 

The  inductive  drop  in  generator  voltage  was  computed 
and  is  also  shown  in  Figure  17.  The  remaining  component  of 
the  voltage  drop  is  that  due  to  plasma  resistance.  By  psing 
this  voltage  drop  together  with  the  generator  current, 
plasma  resistance  and  conductivity  were  determined.  The 
results  of  this  analysis  were  previously  shown  together 
with  the  other  conductivity  measurements  in  Figure  12. 

We  note  that  the  conductivity  estimated  from  the  MHD 
generator  data  is  consistently  higher  than  the  values  deter- 
mined from  upstream  resistance  measurements,  and  close  to 
free-stream  values  obtained  from  eddy  current  measurements. 

It  should  be  noted  that  CG#3  data  may  be  suspect,  since 
fabrication  problems  caused  the  shot  configuration  to  differ 
from  the  calibration  configuration  (the  coil  was  closer  to 
the  magnet,  and  the  magnets  were  slightly  further  apart). 
Generally  the  MHD  generator  data  agree  reasonably  well 
with  the  free-stream  conductivity  values. 
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The  self-excited  generator  output  is  shown  in  Figure 
18.  Voltages  in  excess  of  85  volts  were  reached.  The 
observed  drop  in  the  voltage  at  1 ^ microseconds  is  probably 
an  instrumentation  artifact  since  its  effect  is  not  seen 
in  the  Instantaneous  magnetic  field.  A peak  magnetic  field 
of  O.65  Tesla  was  reached  at  36  microseconds,  an  increase 
by  a factor  of  8.32  over  the  initial  field,  or  18. 3 dB. 

This  is  essentially  the  same  result  as  was  obtained  in 
run  130-3  with  a smaller  diameter  driver. 


The  measured  electrical  conductivity  determined  from 
the  generator  characteristics  was  shown  in  Figure  15  and 
is  noted  to  be  improved  somewhat  over  130-3 . 


3.5  Field  Amplification 

Figure  19  is  a plot  of  the  self-excited  field  amplifi- 
cation of  Shots  130-3  and  130-4.  We  note  from  the  figure 
chat  the  50  percent  longer  time  of  the  flow  of  130-4  did 
not  permit  greater  field  amplification  as  the  trajectory 
in  field-time  space  is  essentially  the  same  as  130-3  with 
the  exception  that  the  absolute  values  of  130-4  are  approx- 
imately 10?  less  than  130-3  at  any  given  time. 


The  conclusion  of  this  plot  is,  that  longer  test  times 
by  themselves  will  not  increase  field  amplification  significantly 
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unless  the  rate  of  field  amplification  for  early  times  in 
the  flow  are  also  improved.  For  example,  if  the  initial 
enhancement  rate  of  130-3  of  1.3  dB/ysec  could  be  maintained 
for  30  microseconds  or  greater,  then  a *10  dB  improvement 
(factor  of  in2)  increase  in  magnetic  field  could  be 
achieved.  For  the  present  conditions  the  field  amplifica- 
tion peaks  out  at  30  to  *10  microseconds  and  the  greater  test 
time  is  not  utilized  to  advantage.  This  demonstrates  that 
future  field  enhancement  pains  must  be  obtained  bv  simulta- 
neously increasing  flow  velocity,  conductivity  and  flow 
duration. 

3.^  Summary  of  Conclusions 

This  new  explosive  plasma  source  can  be  distinguished 
from  most  other  high-performance  explosive  devices,  for 
examnle  the  disc  driver  or  the  Voitenko  compressor,  in  that 
the  geometry  employed  in  the  present  device  permits  a far 
greater  fraction  of  the  pas  tube  processed  initially  by  a 
single  strong  shock.  In  this  way  a higher  average  enthalpy 
plasma  is  obtained  with  the  result  that  one  obtains  a higher 
conductivity  (<j~38  kS/m),  and  the  conductivity  remains  high 
for  a comparatively  long  period  of  time  in  the  flow  (>20ys). 

The  flow  velocity-history  of  this  plasma  source  was 
essentially  identical  when  xenon  was  replaced  by  argon. 
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Bulk  electrical  conductivity  of  the  argon  was  about  70 Jf 
of  the  xenon  values,  and  average  values  from  resistance 
measurements,  about  the  same. 

The  self-excited  MHD  generator  worked  well,  and  per- 
formed in  accordance  with  theory  based  on  measured  field 
coil  and  plasma  parameters.  A maximum  total  gain  of  19  dB 
was  accomplished. 

The  measured  electrical  conductivities  of  the  larger 
source  are  noted  to  be  improved  somewhat  over  the  smaller, 
confirming  that  greater  conductivities  are  possible  using 
the  larger  diameter  plasma  source. 

The  evidence  supports  the  premise  that  higher  velocity 
flow  in  the  channel  can  also  be  obtained  by  the  larger  dia- 
meter source,  but  that  the  higher  reservoir  conditions  may 
also  cause  a ’’pinching  off"  of  the  channel  entrance  and  a 
premature  reduction  of  flow  conditions. 

This  last  effect  may  be  ameliorated  by  simple  geometry 
changes  near  the  reservoir  to  channel  transition  and  will 
be  addressed  in  future  experiments. 
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4 . RECOMMFNDATIONS 

The  progress  and  results  of  the  four  experiments 
designed,  fired  and  analyzed  In  this  period  were  very 
satisfying  and  suggest  additional  fruitful  research. 

We  recommend  a detailed  analysis  of  the  time  history  of 
the  plasma  source  reservoir  conditions  including  the  effect 
of  piston  trajectory  be  made.  We  feel  that  by  understanding 
the  dynamics  of  this  region  we  may  well  be  able  to  tailor 
the  reservoir  dynamics  to  yield  optimum  conditions  in  the 
reservoir  and  hence  downstream  in  the  diagnostic  channel. 


Additionally,  the  premature  closing  of  the  diagnostic 
channel  entrance,  which  was  uncovered  in  the  present  work, 
may  well  have  been  a factor  in  all  of  the  annular  driver 
runs  to  date,  both  smaller  driver  diameter  and  large  driver 
diameter.  Amelioration  of  this  effect  may  produce  significant 
gains  and  should  be  investigated. 


We  recormend  investigating  a plasma  source  loaded  with 
air,  as  the  driver  gas,  as  preliminary  estimates  suggest  that 
the  electrical  conductivity  of  air  at  these  very  high  den- 
sities may  well  be  similar  to  that  measured  using  argon  and 
xenon.  An  experiment  such  as  this  confirming  this  hypothesis 
would  be  of  interest  to  and  have  ramifications  for  other 
scientific  areas  of  interest  to  the  Navy. 
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We  recommend  investigating  different  configurations  of 
conductivity  gages  to  resolve  the  wide  variation  in  measured 
conductivity,  i.e.  to  determine  if  the  variation  is  due  to 
the  variation  of  plasma  properties  with  position  in  the 
channel  or  to  experimental  resolution,  and  a detailed  com- 
parison of  the  conductivity  time  histories  for  clues  as  to 
the  details  of  the  complex  behavior  of  the  flow. 

And  finally,  we  recommend  additional  effort  to  refine 
the  self-excited  generator  diagnostic  section  including 
modeling  of  the  output,  in  order  to  resolve  the  self-excited 
venerator  performance  Unit  problem. 
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